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NONPLANAR WING LOAD-LINE AND SLENDER WING THEORY 


By John DeYoung 

Vought Corporation Hampton Technical Center 

SUMMARY 

Nonplanar load line, slender wing, elliptic wing, and infinite aspect 
ratio limit, loading theories are developed. These are quasi two-dimensional 
theories but satisfy wing boundary conditions at all points along the non- 
planar spanwise extent of the wing. These methods are applicable for 
generalized configurations such as the laterally nonplanar wing, multiple 
nonplanar wings, or wing with multiple winglets of arbitrary shape. Two- 
dimensional theory infers simplicity which is practical when analyzing 
complicated configurations. The lateral spanwise distribution of angle of 
attack can be that due to winglet or control surface deflection, wing twist, 
or induced angles due to multiwings, multiwinglets, ground, walls, jet, or 
fuselage. In a quasi two-dimensional theory the induced angles due to these 
extra conditions are likewise determined for two-dimensional flow. 

Equations are developed for the normal to surface induced velocity due to a 
nonplanar trailing vorticity distribution. Example application of these 
methods are made for a rectangular wing with wingtip winglets. Compared 
with a planar wing, induced drag is reduced 5, 13, 66 percent for aspect ratios 
of 0, 10.9, and infinity. The same results apply to a 30° swept wing with 
0.29 taper ratio and smaller aspect ratio. 

INTRODUCTION 

Theories for multiplane systems and wing endplate theories span the 
history of aerodynamic theory. Early analytical solutions are detailed in 
reference 1. Theoretical development tended to follow the needs of aircraft 
design. In the beginning aircraft were multiplanes including biplanes 
and tri planes. Two-dimensional airfoil tunnel tests in Europe were 



traditionally models of sections with endplates. These required corrections 
from wing with endplates to two-dimensional section data. In reference 2 
by solutions in the Trefftz plane the induced drags are determined of many 
wingtip configurations. These include vertical plates at various span 
stations, circular bodies, and array of fins. Some of these wingtip 
configurations were computed with electromagnetic analogy methods. Winglet 
concepts advanced by Whitcomb are presented in reference 3. These are 
uniquely designed angled wingtip fins for minimizing induced drag and 
improving tip flow. Reference 4 has a chapter on nonplanar lifting-surface 
theory, however, only lift curve-slope data of a V wing and near ground 
data are presented. 

A recent theoretical parametric study of wing winglet configurations 
is made in reference 5. It is concluded in this reference that for the same 
root bending moment, a winglet provides a greater induced efficiency incre- 
ment than does a tip extension. This study provides an abundance of winglet 
data with a wide range of geometric values of wing and winglet. A summary 
statement in reference 5 is 'This study provides sweeping confirmation, for 
a wide range of wings, of the recommendations of Whitcomb in NASA TN D-8260'. 
The winglet concept has been pioneered by Richard T. Whitcomb. 

Present day design objectives are to interrelate aerodynamics with 
structural and weight characteristics and determine possible improvements 
in aerodynamic efficiency. As an example, a planar wing design is 
established. This wing has a certain root bending moment, lift, surface 
area, and induced drag. The designer determines that a wing winglet 
configuration with the same bending moment, lift, and surface area, has 
less induced drag. The designer may choose a highly staggered biplane 
winglet configuration with same as single wing structural weight, lift, 
and surface area, and determines that induced drag is the same, however, 
fuselage bending moment is greatly reduced with the staggered biplane 
configuration. 
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Objectives of the present study are to develop nonplanar wing loading 
methods which are viable for taking into account interference effects 
from other wings or from other vortex systems. Nonplanar means that a 
wing or vortex system extends vertically as well as laterally. 

Because nonplanar wing theory combined with other nonplanar vortex systems can 
become analytically complex and computationally limited, a nonplanar load line 
method will be developed. A planar wing load line theory is developed in 
reference 6. Load line theory is a lifting surface theory confined to a line, 
for arbitrary wing chord distribution, quasi two-dimensional, and solutions 
result in spanwise loading aerodynamics but no chord loading detail. This 
theory satisfies all boundary conditions along the span of the wing. Analytic- 
ally, nonplanar load line principles can be made to extend and expand into 
, nonplanar slender wing and elliptic wing analyses which, in addition, have a 
higher order of direct solution uniqueness. Nonplanar wing theory starts by 
developing generalized loading integral equations which relate the induced 
velocity normal to the surface to the wing loading, or a sum of wing loadings 
if for multiwings. 


SYMBOLS 

A aspect ratio of nonplanar wirtg, 4se^Se 

Ap aspect ratio of planar wing, bp2/S 

ap,an* Fourier coefficient of spanwise loading along the wing, equation 

(15) 

bp wing span of planar wing 

Cg-j induced drag coefficient, D-j/qSe, equation (69) 

Cl lift coefficient, L/qSg, equation (54) 


3 



Cl 

Cmbtj 

Cmbr 

c 

Cav 

Cr)j Cn* 

Cp 

Can 

Di 

Cnn* 

e 

G 

^nn* 

^n* 


wing rolling moment coefficient, rolling moment/qSe2se» 
equation (61) 

wing running bending moment coefficient about lateral point njj, 
Mbb/qSe2se, equation (56) 

wing root bending moment coefficient, Mbn/qSe2Sg 
wing chord 

wing average chord, Se/2Sg 

equals Cp with n or n* substituted for p 

wing chord Fourier coefficients, equation (35) 

Fourier coefficients for product of wing chord and angle of 
attack conditions, equation (35) 

wing induced drag 

nonplanar induced drag influence coefficients, zero for planar 
wing, equation (67) 

induced drag efficiency factor, equation (116) 

dimensionless circulation, r/2sgV, equation (14) 

wing chord influence coefficients, equation (35) 

nonplanar wing chord influence coefficients, zero for planar 
wing, equation (35) 


4 



surface loading factor, equation (24) 


wing lift, also wing rolling moment, equation (61) 

nonplanar wing lateral -vertical plane influence function, 
equation (32) 

Li integration nonplanar wing spanwise loading influence 
function, equation (33) 

odd integer, number of terms in quadrature formula, equation 
(18) 

wing bending moment about a specified span station 

integer denoting number of unknowns and equations, equation (34) 

integers of spanwise loading Fourier coefficients and influence 
coefficients, odd only for symmetric loading, equations (15), 
(34), and (35) 

dynamic pressure, ^ pV^ 

surface area of nonplanar wing, includes wing area in yz-plane 
surface area of planar wing 

spanwise coordinate along nonplanar wing surface, figure 1(a) 

total semi span distance spanwise along wing surface from wing 
root to wing tip, figure 1(a) 


free stream velocity 



w 


w 


y 

z 


a 


a 


C 


r 

Y 


e 

6 

n 

'icp 

0, 0Q 

X 


induced velocity normal to wing surface or trailing vorticity 
sheet, equations (5) or (6) 

integers which denote span stations, equation (41) 

lateral position of spanwise wing station s 

vertical position of spanwise wing station s 

aircraft angle of attack 

wing geometry angle, equation (28) 

induced angle normal to wing surface, equation (27) 

wing circulation 

nonplanar wing dihedral cant angle, measured from y-axis, 
figure 1 (a) 

wing twist angle 

dimensionless vertical coordinate, z/Sq 
dimensionless lateral coordinate, y/Sg 

spanwise station at which wing running bending moment is taken 
lateral semi spanwise center of pressure location 
angles used in derivation of induced velocity, figure 1(a) 
wing taper ratio, tip chord/ root chord 
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planar wing winglet taper ratios, equation (75) 


Awo ratio of winglet root chord to wing root chord, equation (75) 

A sweep angle of wing half chord line 

K integers for quadrature formula, equation (18) 

p dens i ty 

<P spanwise trigonometric coordinate, cos'^a, equation (15) 

a dimensionless spanwise coordinate, s/s 0 , equation (14) 

ip winglet toein angle, positive for toein deflection, equation 

(28) 

Subscripts: 
av average 

b bending moment, also bending moment about lateral point 

cp center of pressure 

e nonplanar wing tip 

i, j ith and jth vorticity sheets, equation ( 6 ) 

n,n*,p integers, equations (34) and (35) 

0 spanwise station where winglet begins, that is junction of wing 

and winglet 
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planar wing 


wing root 
vorticity sheet 

normal to surface induced velocity or downwash point 
integers, equations (41) and (71) 


integers, equations (18), (21), and (22) 



DEVELOPMENT OF THEORY 


Nonplanar loading theory differs from planar theory by the difference in 
velocity induced normal to the wing surfaces. In load line theory this 
induced velocity is simply that due to the trailing vorticity system. Using 
the nonplanar induced velocity combined with the planar load-line theory 
of reference 6 will provide a nonplanar load-line theory. 


Induced Velocity Integral Equation of Nonplanar Vorticity Sheet 

Single vorticity sheet . - A single nonplanar trailing sheet extending 
longitudinally from zero to infinity is shown in figure 1 (a). The normal 
induced velocity at point (y^^, z^^), by Biot-Savart law, is given by 

— dr 

dw(yw.zw) = cos 0 

> 

where r = [(yy - + (zy - ^ 

Referring to figure 1 (a), the following relations can be developed: 


6+'2-+Yw“f‘'''9o> 0“So~Yw 

cos Yw “ COS 00 COS Yw sin 0 q sin yw 

cos 00 = (yy-yw)/>". sin 0q = (zy-z„)/r 

tan Yw = Yw = [1 + YW = (^)^^ [1 + 

Let s be distance along the surface of the vorticity sheet measured from 
midpoint of the sheet. Then the differential relation between s and y is 

ds = [1 + (0)2]1/2 dy 


h (2) 


( 3 ) 
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With the 9, 00 * and yw values of equation (2) inserted into equation (1), the 
incremental downwash at point yy^, Zy^ due to a nonplanar trailing vorticity 
sheet is 


dw(y„,z„) = - [yv - yw + (^v - Zw) tan ^ ds 


[1 ^ (f)y 

Integration of equation (4) leads to 


1 r, . ,dZx21-l/2 yy~yW + (Zy-Zy^) (jy)yy ^ 

(ay)w3 ^ ds 


(4) 


ds 


w 


fv 7 ^ = _ I_ n + yV~yw’*'(^V“Zyy)(j j 

(yw.zwj - 4TT Ll + 7 ^2-— ^ d? 

7v-yw) + (zv-zw) 


ds 


(5) 


where the integration is taken along s, that is a line integral along the 
vorticity sheet. The integrand variable y^ is related with s by the expres- 
sion given in equation (3). The point at which normal velocity w is 
determined is defined by and z^. The nonplanar vertical displacement of 
the wing given by z is an arbitrary function of y. The derivative (dz/dy)y^ 
is the slope at z at the point z^, and is constant in the integrand of 
equation (5) as are also y^ and z^. 

Multiple vorticity sheets . - An example of a pair of vorticity sheets is 
shown in figure 1(b). The total induced normal velocity at sheet j is not 
only that due to sheet j but also that due to sheet i, or for many sheets with 
i = 1, 2, . . . In equation (5) the induced velocity normal to the slope 
line (dz/dy)y^ is determined at point (yyy, z^) which is a point on the 
vorticity sheet represented by the function Zy{y). However, the point 
(yw> Zyy) is arbitrary and can be on any curve z(y), that is can be different 
than the Zy(y) curve. Let Zy\ represent an N-group of vorticity sheet 
curves, then the normal velocity along the Zwj curve obtained by reapplication 
of equation (5) is given by 


wj(ywj.Zy^j) - n + 


N 

%r 

1=1 -Sei 


yvi-ywj-^Zvi-zwj)(^)wj ^ ds„,l(6) 
(yvi-ywj)" (Zvi-Zwj)^ 
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where j = 1 through N. With equation (6) the induced normal velocity can be 
determined along the span of N abitrarily shaped vorticity sheets or wings. 
Some examples of various shapes that can be analyzed are shown as follows: 




o — c 



Induced Velocity of Wing with Winglet 
For a single wing with winglet the vorticity sheet appears as follows: 



For this nonplanar geometry the coordinate relations are 


in the range 

-yo^y 6 yo 


in the range 

yo ^ y ^ye 


in the range 

-ye y ^ -yo 


z = 0 



y = s 


z = (y-yo)tan Y 
0 = tan Y 

y = So + (s-so)cos Y 
z = (s - So)sin Y 


z = (-y-yo)tan y 

(7) 

^ = -tan Y 

y = -So+(s+So)cosY 
z = -(s+So)sin Y 


With equation (7) values inserted into equation (5) the velocity normal to 
the surface is given by 
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for 0 ^ Sw ^ s 


w(sw) = 


'kf 

-So 


Sn ^ 
dsv 


ds\ 


$v - Sw 


V 

-So 


idr 


^ ^ ®[SvCOSY -Sw+So(l-COSY)]^ydSv 

[SyCOSY -Sv;+So(1-COSy)]^ +(Sy-So)^Sin2Y 


jj^^O[Sy COSY -Sw-So(I-COSy)] dSy 


( 8 ) 


4ti^ 

-S( 


[sy COSY -Sw-So(I-cosy)]^ ( Sy+SQ ) ^ si’ii^y 


and for Sq ^ ^ Se 

w(Sw) = 


s dr 

1 p ° [SyCOSY-Svv+So(1-COSY)]^^ dSy 

[sy-SwCosY -So(I-cosy)]^ +(Sw-So)^sin 2 Y 


-hf 

-Sn 


TY~ 

in^Y ^ 
So 


-e ^ ds 
dsw 


Sy - Sw 


[-(1-2cOS^Y)Sy-Sw-2So(1-COSY ) COSY ] dSy 
Vg [(sy-Sw)cosY -2 sq (1-cosy )]^ + (sy+Sw)^ sin^y 


(9) 


The integral with the (sy - Sy^)"^ term in the integrand contributes a 
large part of the induced velocity, because it results from the vorticity 
nearest the point at which the velocity is determined. This integral is 
analytically simple to obtain if the integration is taken over the whole 
wing; that is, from -Sg to Sg. Equations (8) and (9) give the same result if 
the (Sy - Syy)"^ integral is added for extending the integration and subtracted 
from the more complicated integrand. With this adding and subtracting and in 
dimensionless notation, equations (8) and (9) simplify to 


for 0 ^ Qq 


“i (%). 




(10) 
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( 11 ) 


and for oq ^ % — 1 




where 


L± = - IT 




(oy ± 0o)(ay + aw ± 20n)(l - cos y) 

)-^ + 2(oy ± ao)(aw ± Oo)( 1^0S'Y)J (CTy-Oy,) 


( 12 ) 


I (1-2 cos^Y)qy + ow •*• (1~9Q^ - 1 

0 C(ay-cTy^)cos Y - 20 q (I-COSy)P + (cry+ay,)^sin^Y cry “ 


(13) 


The dimensionless values in equations (10) through (13) are defined as 
f ol 1 ows : 


a = 


Se 




= ^ 


G = 


2SgV 


w 

Oli = TT 


(14) 


With equations (10) and (11) the induced angle normal to the nonplanar 
wing surface can be determined for an arbitrary loading G. This is for the 
nonplanar wing or wing-winglet defined in equation (7). 

Integration with G represented by a Fourier series . - Let the wing 
loading coefficient along the s or a coordinate of the wing be given by 

N N 

S(4>) = ^ an sin n4> , ^n an cos ni> (15) 

n=l n=l 


where a = cos<f>. Then the first integral in equations (10) and (11) becomes 

1 A dG . . N IT 

_ 1 r day ^ l_r d»v ^1 Vn a cos n4»y dj»w 

Oy - Oy, ^ COS(J)y - COS^ 2iir 2_ '\y COS(})y - COS(j)y^ 

-1 0 n=l 0 


Since 


IT 



0 


COS njiy d(|>y 

cos (|)y - COS (ji^ 


trsin nj>w 
sin 


(16) 
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then 


1 dG . 

l_r d£vL^= 1 

2-nJ Oy - CTy 2 sin 


N 

^n an sin ni})^ 
n=l 


07 ) 


When G is expanded in a Fourier series the L integrals become complicated; 
however accurate integrations can be obtained by using a simple quadrature 
formula given by 

^2 M 

f' f(*)d<,= [f(n) * fM H. (,8) 

P=1 


where = <})i + ^ M is an odd integer. Then the L_ integral 

of equation (10) becomes 



^0 


L 


i’o 

^d0y = ^ L_(<|.v,<^w)^ d<()v 
0 



L_ { <l>y s <f'w)cos n(f)y d^y 


N Ml 

■ ( 0 , (jiw )+L_ ( ^ 0 » '/’w ) cos n<|>o+2 ^ L_(<})ip,(j),^) 

n=l ^ p=l 


cos n<j)i^JJ* (19) 

where (Jiiy = and where examination of equation (12) shows that 

L_ (<l>o»<f>w) ^ Similar integrations are made for the other L integrals 
which appear in equations (10) and (11). With the equation (19) procedure 
and with equation (17), equations (10) and (11) for symmetrical loading 
(n = odd) reduce to 

N 

«l(*w) = Z " "n (riT^^ + L„„) (20) 

n=l 

odd 
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where for 0 ^ ^ ctq or 4 >q^ (Jv ^ ^ 



( 21 ) 


where c^iy = <p2u = ^ - 


■H.to_ 

M2+I 


and for oq ^ 1 or 0 ^ ^ (j>Q 

‘-m = r-l- 4 i'-»g.tw)co . s nfa + J l.(+3v.+w)cos n*3„l 


Ms 

I 

y=l 


M 4 


■*■ 2 tt (M4+I ) o(^»<t>u) 4.^ Lq(,J)i^^ ,<j)^^)cos n(J)4jjJ 

y=l 


( 22 ) 


where <}> 3 y = <|>o + . <|) 4 y = ^ - and Mi, M 2 , M 3 , and M 4 are 

odd integers. 

In equations (21) and (22), the L+(0^,<(.w) and Lo(<^u,c(.w) functions are given 
by equations (12) and (13), in which = cosi})y, aw * cos and Cg ■ cos (>o* 


Load-Line Theory 

Load-i.ine theory is developed in reference 6 for the planar wing. In 
this theory the boundary conditions are satisfied at all points along the 
span and an aspect ratio factor is introduced, however, the theoretical 
simplicity of the lifting-line method is maintained. In the present develop- 
ment, this load-line theory is extended into a load-line theory for nonplanar 
wings. 

The integral loading equations start with the relation that section 
loading is given by 
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pVr = 2ir - ka-j)qc 


( 23 ) 


or 



kai ) 




where is the incidence angle between free stream velocity and wing chord, 
a-j is the induced angle due to wing wake vorticity sheet or multiple 
vorticity sheets extending longitudinally from zero to infinity, and k is 
a lifting surface downwash factor developed in reference 6. The factor k 
originated in the study presented in reference 7, where it had the value 
(A + 4)/ (A + 2). Later work showed that a more accurate value for k is 


_ A + 3.79 
A + 1.895 


(24) 


Using the dimensionless definitions of equations (14) and (15) and defining 
nonplanar aspect ratio and average wing chord as 

A = , Cav = If-g* Se = surface area of nonplanar wing 

2se = wing tip to tip perimeter length 

then equation (23) becomes 

N 

^ an sin n<j> = f ^ (“c " l<“i) (26) 

n^ ^ 



Equation (26) is the form of the equation for lifting-line theory 
solutions except here the k factor is introduced. In lifting-line theory 
k=l . The planform geometry terms c/cav and are functions of (}>, and d-j = 
a-j((j>, n, ap). By specifying (j> at m spanwise stations, .m linear simultaneous 
equations result with m unknown ap's. These are solved for ap. These 
values of ap represent the solution for which the boundary conditions are 
satisfied at a finite m spanwise stations. 
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In load-line theory the boundary conditions are satisfied at all points 
along the wing span. This is possible by a spanwise integration of the wing 
chord, angle of attack, and induced angle. In equation (26) ai can be 
represented by the function 


«i = ai ((}>, n*, ap*) 


Then the expression on the right side of equation (26) is only a function 
of 4> and n* but not of n. Then by Fourier theory the Fourier coefficients 
are given by 


an = 



- ka^*((j), n*, an*)] > sin n(|) d(}> 


or 


an = 



^ — ac Sin n(j) 
^av 



((j),n*,an*) sin n(j) d<() 


(27) 


for n-1,2,3, . . • 


For a single wing, ai in equation (27) is given by the ratio of equation 
(5) to V. For multiple wings, ac in equation (27) is added to equation (6) minus 
equation (5). In addition, all the terms in equation (27) become those for 
the jth wing; that is, A j , (c/cav)j. ^c j , k j , and ajp, which are the Fourier 
loading coefficients for the jth wing. The wing geometry angle ac is made up 
of the effect of wing or aircraft angle of attack a on the nonplanar wing, 
the wing twist e, and the toein-toeout angle i|) of the nonplanar wing. Wing 
twist and toein angle with the wing uninclined are defined as angles between 
the wing chord and free stream velocity, positive in the direction for loading 
increasing on the upper wing surface. For the nonplanar wing the effective 
angle of attack at a wing section is a cos y. Then 


oiQ = CL COS Y e ^ 

where tan y = (^)^, cos y= [1 + (^)J 
and e = e((j)) , \i> = ij((())) 




(28) 
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Thus the solution of the nonplanar wing loading integral equation reduces to the 
load-line solution for determining the an Fourier loading coefficients of 
equation (27). It is simply a wing geometry integration which includes the 
wing chord distribution along 0 or (j> and the angle oq distribution given in 
equation (28). The second integral involves the integration of w (y^, Zw) " 
Vai((i)v^) of equation (5). In a dimensionless equation (5) the nondimens ional 
derivative of circulation becomes an n* series given by equation (15) but 
with n = n*. and <p = ipy. 

Load-line method for single wing . - The induced angle in equation (27) 
can be expressed as that due to a planar wing plus that due to the difference 
between nonplanar and planar wings, as was done in equations (10) and (11). 

Let a-jp represent the induced angle due to a planar wing and aio that due 
the effects of the nonplanar condition. Then «i = aip + aio- Then, in 
nondimensional terms, equation (5) can be written as 


(<(>w) = “ip (<f>w) “io (4>w) 


(29) 


where with equation (17) 
“ip ( 4’w ) ~ 



N 

n*an* sin n* ((>m 
‘ 2 sin <j>w 


and with equation (15) 

“io (‘f’w) 




doy 


N 

y ’^*^n*^n*w 

n*=l 


(30) 


(31) 


where 


nv ~ nw + (?v “ ^w)(^)w 1 

+ (Cv-C„)^] ■ "V - Ow 

which approaches zero as ny n^,, 


(32) 


18 



1 


Ln*w = 

0 

n = y/Se» C = z/se, o = s/se = cos (|), and n. C, and a by 
equation (3) are related by dcr = [1 + dn 


(33) 


With equations (29), (30), and (31) inserted into equation (27), the loading 
integral equation for the single nonplanar wing becomes 

N N 

^ n*Inn* ®n* " ^ ^ '^*l^nn* ®n* (34) 

n*=l n*=l 


_ ir irk 

«n " A “ 2A” 


n = 1 , 2, 3, 


N 


where 


i TT 

Can = I/' «c sin n<). d<|> 


^nn* “ - 


c sin n’ 


i sin n<t> d4> = 


ir ] Cgv Sin <|) 
o 


2 r C P^'' 4 

^ p<^ d.1. = - 

7T 

•^n* = I ^ l-n*w sin M d<j. 


n+n*-l 

^P 

p=|n-n*l + 1 


(35) 


where the I^n* - Cp relation is derived in reference 6. The first two 
terms on the right side of equation (34) are the same as. those for the 
planar wing since then K^n* = 0 because Li = 0. The solution for ap involves 
N equations and W unknowns of ap or an* which are solved as linear simultaneous 
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equations. The coefficients Cp and Cap are integrations of the boundary 
conditions, that is, the wing chord distribution and wing angle of attack 
distribution along the spanwise (J> = cos"^ a parameter. The cp's are the 
Fourier coefficients of the wing planform which wing in a Fourier series 
function can be represented by 

oo 

= y Cp sin p(|) (36) 

^av L- 
p=l 


Similarly Can's are the Fourier coefficients for the product of wing chord 
and geometric angle ac. The angle which the wing chord makes with free 
stream velocity is given in equation (28) where a is constant with respect 
to 4> or a. Equation (34) applies to a single nonplanar wing of any vertical- 
lateral shape since ^ = ?(n) or z = z(y) are arbitrary functions. In addition, 
the wing chord distribution with unswept midchord line can be for any plan- 
form shape, and the wing twist or winglet angle distribution are arbitrary 
functions. 


Symmetric spanwise loading . - This loading results from a spanwise 
symmetric wing planform with spanwise symmetric distribution of angle of 
attack across the span. With these conditions only odd numbered values of p, 
n, and n* of Cp, Cap, and Knn* have values in equation (35). For symmetric 
loading solutions it is convenient to evaluate the ai coefficient in terms 
of ratios of ap/ai. For n = 1 equation (34) becomes 


an 


I. Tfk T 

A ^al - 2A 


ai 



N 



n*=3 

odd 


n*Iin*an* 


N 

^n*Kin*an* 

n*=3 

odd 


(37) 


4 

Since In ~ Im* " ^n*’ ^ni ” ^n 


then , _ ■^^ai 

^ 

A + 2k [1 + jKii + J ^ (cn* + 2Kin*)n* (38) 

n*=3 

odd 
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Dividing equation (34) by equation (38) gives 

N 

ai 


= [n|li(l+5Kii)] - f^(cn+2Kni)+ J^2[|^(cn*+2Kin*)-Inn*-2Knn*]n* 

(39) 


n*=3 

odd 


for n=3,5,7, ...N 


where the coefficients Cotns Inn*> l^nn** (note that p can be substituted 

by n or n* in the Cp integral) are given in equation (35). Only the 
coefficients for odd numbered values of n, n*, and p are needed for the 
symmetric loading solution. Equation (39) has (N-l)/2 equations in n with 
(N-l)/2 unknowns of an/ai or an*/ai which are evaluated by solving as 
simultaneous linear equations. The Cn and Cn* coefficients are the Cp values 
given in equation (33) but with n or n* substituted for p. 

Planar-wing winqlet solution for symmetric spanwise loading . - The Fourier 
loading coefficients a^ are obtained from the solution of equations (39) and 
(38). For this nonplanar configuration the function of equation (32) is 
formulated as and Lq functions given in equations (12) and (13). The 
numerical integration of these functions results in the Lpw expressions given 
in equations (21) and (22). From equation (35) the Knn* coefficients are 

T7 

Knn* = I ^ Ln*w sin n<t> d<j> (40) 

0 


where for this nonplanar configuration the are given by Lnw of equations 
(21) and (22) but with the substitution of n* for n. Equation (40) can be 
integrated numerically by using the quadrature formula of equation (18). Let 
i> = <t>w ~ wtt/(M+1), then 

M 

Knn* = fIt I Ln*w sin n<^wl 

w=l 
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where M is an odd integer. Since the chord distribution spanwise is 
symmetric and for odd n sin n(j>w is symmetric, and the induced angle 
contribution of Lp*y; is symmetric, then 

M--1 

2 c c 

Knn* = t(^)M+l '-"* M+1 sin n f + 2 ^ (— )^ Lp*w sin n<j.w] (41) 

~ 2 ~ ~ 2 ~ w=l 

where the (M+l)/2 subscript denotes midwing conditions, that is at = tt/ 2 
or = cos = 0. The chord term (c/cav)w ’’S that at span station ^ = 
wtt/(M+ 1). Lp*w is with (j)^ in the <|)w range of equations (21) or (22), in 
which n is changed to n*. In equations (12) and (13) = cos [wir/(M+l)]. 

To minimize numerical integration error c|)q should be chosen between the 
(fiw points, or (j>o ^ wir/(M+l). 

The planar-wing winglet ap coefficients for symmetric spanwise loading 
are obtained by the solution of equation (39) with the Kpp* of equation (41) 
and the Cap, Inn*> Cp of odd n, n*, p integers in equation (35). Example 
values of Lp*^ with (j>o = 5ir/32 radians have been computed using equations 
(21) and (22) and are presented in table I. These values are for M = 15 
with (fiw = wir/16. 


Slender Wing Theory and Elliptic Wing Theory 
Elliptic wing and A-»«. limit solution . - For the elliptic wing 

c = Cp (1-CT^)^^^ = c^ sin (j), $e = ^ Cr^e* ~ ^ 


_Se c _4/n_?xl/2_4_..„ 

- ?f- - 2T cr. ^ - 7 (1-^^ ) " 7 


(42) 


With this c/cav equations (26), (29), (30), and (31) are combined to give 
N N 

^ (A + 2 nk) ap sin n <}> = 4 ac sin<fi - 4k ^ n* ap* Lp*^ sin <f> (43) 


n=l 


n*=l 


22 



then the Fourier loading coefficients for the elliptic wing planform are 

N 

= A+Iiik ^ Y. 

n*=l 


for n ■ 1 , 2, . . . N 

where Cap' and Kpp* are the xoefficients given in equation (35) but with c/cav * 
(4 /it) sin (j), that is 


Can 



sin <f> sin nc{> dc|). 


Knn* 



sin (p 


sin 


nc{> d(j> 


(45) 


For symmetric loading on the elliptic planform, equation (44) becomes 


TTC 


ai = 


al 


A + 2k(l + ^ Kii + ^ ^ n*Kjp* 1^) 


n*=3 

odd 


S (46) 


N 

Kn-Kni)+^k y (^i^*-Kpp*)n* 

£~ v-ai 

n*=3 

for n = 3, 5, . . . N odd 

As aspect ratio approaches infinity, the solutions for ap reduce to direct 
relations between ap and the wing geometry coefficient Cap. For A -> «> 
equations (38) and (39) become 

TTCaj 


ajl Cqp 
ai Cqi 

where Cqn is given in equation (35), or for the elliptic wing by equation (45). 


(47) 


Inl] 

ai J 


an = 1 

3i A+2nk 




A ^ + 

^al 


Can + iion 

' 'TT r' r* 


' 7 T C 


al 


Cgn 
Cgi 
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Slender wing . - When aspect ratio approaches zero, then^ for any wing 
planform with trailing edge cutout aft of all points of swept leading edge, 
as shown in reference 8, equation (44) for A 0, k 2, becomes 

N 

®n ~ ^ (^an ^ n* Kpn*ap*), n = 1, 2, . . . N (48) 

n*=l 

where Can and Knn* are given in equation (45). With symmetric loading 
equation (46) for a A 0 wing is 


ai 


1 

N 

4 + 2ttKii + 2 it ^ n* K 2 n* 

n^=3 

odd 





(49) 


for n * 3, 5, ... N 

Another derivation of nonplanar slender wing theory is to note from 
equation (26) that the boundary conditions are simply that ac = ka^ = 2ai 
as A ->• 0. From equations (29), (30), and (31) 

N N 

“i " 1 . sin" J " ®n sin n<^ + ^ n* a„* Ln*^ (50) 

n=l n*=l 


then, with the condition ct^ » 2 ol ^, equation (50) can be written as 

N N 

Y n aj 2 sin n cf> = ac sincfi - 2 ^ n* a^* sin 

n=l n*=l 
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from which the Fourier coefficients are 

2 2^ 
naf, = — / oig sin <j) sin n<() d<j) - 2 Y n*an*(;|j^Ln*^^ sinc|) sin n(|) dij)) 

a. T -^0 

0 n*=l 

which in terms of Co,n and Kp^,* of equation (45), becomes 

N 

^n ~ 4 ff (^®n ” ^ ^ ®n*^ ,n-l,2,...N (51) 

n*=l 


which is the same as that of equation (48). 


Force and Moment Coefficients 


Spanwise loading and lift coefficient . - With the Fourier loading 
coefficients ap determined by equation (15), the spanwise loading 
distribution along the spanwise a or coordinate is 


G = 


r 


N 

^ ap sin n<() 
n=l 


(52) 


where <p = cos~^a, and the loading is normal to the wing surface in the 
yz - plane. 

The lift coefficient is given by 


Cl = 



siniji d<() 


(53) 


where is the vertical component of the loading coefficient G which 
itself is normal to the nonplanar wing surface". G^ is given by 

Gz = 6 cos Y 
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where as in equation (28) 


cos Y = [1 + 


then 


ir N IT 

C|_ = ~ ^ J ^ <fi d(f> = A ^ a^ f cos Y sin({) sin n(|> d(j) 

® X n=l '^0 


ttA 


9l 


N ^ 

- A ^ ap j ( 1 - 
n=l ^0 


cos y) sin sin nij> d<J> 


(54) 


For a planar-wing winglet, y is constant between 0^ cji ^ (j)o» 

and zero between <t>o ^ ir-<}>o» then equation (54)for symmetric loading 

reduces to 


Cl = - A(l-cos y) 



lin^) ^ 


Ia„[ 


sin(n-l )(fio 
n-1 


n=3 

odd 


sin(n + 1)^0 11 

n + 1 Jj 


(55) 


Wing Bending Moment . - A wing running bending moment along the wing 
span coordinate a or ^ can be formulated for a nonplanar wing by use of the 
following geometry: 
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The incremental bending moment at point (yb>Zb) is the incremental loading 
pVrds times the perpendicular moment arm between points (yt, zjj) and P. 
That is 


dMbb = pVr [(y-yb)cos y + {z - zb)sin yl ds 


With integration and in dimensionless terms, the wing bending moment 
coefficient at point (nb»?b) 

Mbb A 

^bb = - qs ^ gse " 7 / y + (?-?b)siny]G da 


^b 


A / ^ 

= ^ / [(n-nb)cos Y +(c-<;b)s1n y]G sin 

0 

% J 


[(n-nb)cos Y + (?-Cb)sinY]sin(() sin n<|) dcf. 


where G is given in equation (52) and the coordinate relations and y are 
defined in equations (28) and (33). 

For a planar-wing winglet the coordinate relations with nb ^ no af’e 
given by 


(56) 


n-rib = 
?-?b = 0 
Y = 0 


0 ^ a o, 


^0 • 


n-nb = ao-a| 3 +(a-ao)cos y 
C-Cb = (a-ao)sin Y 
Y = constant 


Lao^cr^l (57) 


then equation (56) simplifies to 
A 

^bb ~ 1 / (cr-0|3)Gda ^ ^ i C *^“®0 (0o~^b)coSY]G da 

<^b *^0 


27 


II nil II 



where y is constant in the integrand. For symmetric loading with G in the 
Fourier series expansion, for the planar-wing winglet and <f)fj - (|)os this 
integrates into 


Cmbb " 8 J 


an 


n=l 

odd 


j sin(n-2)4ih sin(n+2)(j)h 
1 n-2 " n+2 


-2 cos (|)t) 


■^b-sincf>bCOS<|)b, n=l 

sin(n-l )(i)h _ sin(n+l )c|)b 
n-1 ~ n+1 



+ 2(l-cos y)(cos (j>b-cos (|)q) 


■<j>o“Sin (j)o cos (j)o, n = 1 


sin(n-l )4>n sin(n+l )j>o 


n-1 


n+1 


, n >1 


( (58) 


where Cmbb is the wing bending moment coefficient at span station (|)b due to 
loading outboard of ({)b. The wing root bending moment is with i()b = J, then 
equation (58) reduces to 


^mbr “2 ^ | 


n=l 

odd 


1^2 _4 "" “^(l— COS y)cOS <J)q 


(j)o-sin (JjqCos (f>o, n=l 


sin(n-l) (|>n sin(n+l )(j>n _ 

n-1 ■ n+1 ’ " ' 

(59) ' 

The coordinate relations for running bending moment of the planar-wing 
winglet at points on the winglet, nb — no> 


n-nb = (o - cTb) cos y 

e-Cb = (p - Pb) sin y 

y = constant 


then using equation (56) as before 
1 

Cmbb ^ (p-Pb)C da 

Pb 

N r 

_ A r J sin(n-2)d,b sin (n+2)<j>b o ,, 
8 / 1 n-2 h+2 ' ^ 

n=l I 

odd 


>b-sin<()bCOS(j)b, n=l 

sin(n-l )(j)h _ sin (n+1 )j)h 
n-1 ~ n+1 

n>l 

(60) 


28 



where the bending moment is at span station <j>b ^ <|>o> that is, ({ib is on the 
winglet. 


Wing rolling moment . - Rolling moment is the same as the wing root bending 
moment except the spanwise loading is antisymmetric and the integration spans 
the wing. Thus for rolling, equation (56) becomes 


Cl 


L 

qS02s0 


A 

7 


0 


TT 

(n cos Y + ? sin y)G sincj) d(j> 


(61) 


where here L denotes rolling moment and Ci is the rolling moment coefficient. 
For a symmetric planform and antisymmetric loading 


Cl = 


N-1 tt/2 

A y a„ r 

n=2 4 

even 


(ncos Y + 5 sin y) sin <j> sin ncfi d(j) 


(62) 


The spanwise symmetric planar-wing winglet specified in equation (57) has 

r 

C] = A / aG da + A j {o-oq + ao cos y)C dt() 
o ao 


then with antisymmetric loading 

C] = ^ a 2 - ^(1-cos y)cos ij)o 


N-1 


I [ 


sin(n-l )(^a _ sin(n+l )(|)o 
n-1 ~ n+1 


(63) 


even 

Induced drag coefficient . - Induced drag is the spanwise integration of 
the product of spanwise loading normal to wing surface in the y, z plane and 
wake normal -to-surface induced angle. Thus 


Di 



Taj ds 


(64) 


29 


I I 1 1 1 II II I 



In dimensionless terms 


1 


Ga-j da = A 

With equations (52) and (50) inserted into equation (65), the drag 
coefficient becomes 

N N N 


G(<j») a-|(<j)) sin<() d(|) 



C[)^. - ^ ^ na^ + A ^ an ^ n*an* Dnp^ 


n=l 


n=l n*=l 


7T 

- f Lp*v^ sincf) sin ncf) dif) 


(65) 


( 66 ) 


(67) 


where Ln*w is given in equations (33) and (32). This Dpn* integral is of the 
same form as the equation (35) Knn* integral and is identical when (2/ir)(c/cav) 
equals sin <j>. Comparing equations (67) and (45), it can be noticed that 
the Dnn* coefficients have an identity with the Knn* coefficients of the 
elliptic wing; that is 

Dnn* = ^ '^n*eniptic " 8“ 

wing ^ 

Examination of equations (33) and (32) shows that Dnn* depends only on the 
wing spanwise vertical displacement, the ? = ^(n) function, and is 
independent of all other planform geometry such as aspect ratio, chord 
distribution, sweep, and wing angle variation. Thus once Dpn* coefficients 
are tabulated for a given ? curve, they are valid for any planform shape. 

For symmetric spanwise loading, Ln*w is symmetric, and also only odd n and n* 
apply; then equations (66) and (67) simplify to 
N N N 

Coj = ^ ^nan+A^ap ^ n*Sn*Dnn* ^^D) 

n=l n=l n*=l 

odd odd odd 
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Ln*w sin (p sin n4> d(|>, for odd n and n* 


(70) 


Dnn* 


u/2 

= !/' 


The Dnn* coefficients of the planar-wing winglet configuration with 
symmetric spanwise loading are determined in the same manner as done for the 
Kpn* coefficients of equation (40). After employing this similarity relation, 
which is simply the substitution of sin (j>w for (2/Tr)(c/cav) into equation 
(40), equation (70) becomes .j 

•^nn* = (Ln* m±V sin n J + 2 V L„*^^sin «j.„ sin n<|.„) (71) 

~t~ w=l 


where cfy, = wir/(M+l), and M, n, and n* are odd integers. For the planar-wing 
winglet Ln*w's are determined from the appropriate (jy, range of equations (21) 
or (22) (in which n is changed to n*). 

Example computations have been made for one planar-wing winglet 
configuration; that is, for <f>Q = 5ir/32 radian or oo - .8819 which is the 
span station at which the winglet begins. For M = 15, then (!>« = wir/16, the 

computed from equations (21) and (22) are presented in table I. With 
these Ln*w values, the Opp* coefficients are determined from equation (71) 
using M = 15. Dpn* coefficients are presented in table II. When 
these Dpp* values are inserted into equation (69), the induced drag coefficient 
depends only on the spanwise circulation distribution or laterally 
surface spanwise loading distribution normal to the surface, and is 
independent of planform shape including winglet shape or sweep. This 
planform is nonplanar and is for a planar-wing winglet with wing tip winglets 
starting at oq = .8819 span station with winglet at a given y angle. 
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PLANAR-WING WINGLET LOADING SOLUTION 
FOR RECTANGULAR WING 

The sample application of the theory developed in the previous chapter 
demonstrates the usage of the method. The example wing configuration is a 
symmetric wing with winglet at each wing tip, the combination being described as 

a planar-wing winglet. Particularly, the rectangular wing and the slender wing 
is investigated in detail. An objective is to show example solutions of 
equations (39) and (49) for the prediction of an/ai Fourier loading coefficients. 

Chord Distribution for Tapered Planar-Wing Winglet 
The chord distribution is 


c = cr - (cr-Co) 0 s a < Oq 

(72) 

= C\fiQ - (cwo"^wtHyE^) ’ oq < a < 1 

where Cp is wing root chord, Cq is wing chord at span station oq , Cy^Q 
winglet chord at span station ao> and c^t is winglet chord at winglet tip 
(ct = 1). The chord distribution of equation (72) leads to the following wing 
geometry relations 


Se 


Cr^et (l'*’^o )*^o'*’^wo (i “'^o > 




(73) 


Cr = 2 

Cgv (1 +Ao)cTo'*'Xwo(1'^^w) (1“Oo) 


(74) 


_9_ = £r_ (1 

^av Cgv 


^ a) » 0 < a < oq 


= C2i_ 
Cav 




^0 


1 a s. 1 


( 75 ) 
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where = Cq/c^ is planar-wing taper ratio, x^ = %t/Cwo is winglet taper 
ratio, and x^q = c^q/c^ is ratio of winglet root chord to wing root chord, 
and a = cos <(), oq = cos <J)q. 


Planform Fourier Coefficients 

From equation (35), the Cp coefficients are obtained from the Integration 


n 



4 

ir 


(76) 


where, for the tapered planar-wing winglet, c/cav is given by equation (75). 
For a rectangular wing with rectangular winglet, equations (74) and (75) 
reduce to 


c ^ 1 

Cav ^wo'^Ci'^wo^ ^®s <j>o 


(}>o 2 <t> s 


TT 

2 


, 

^wo'*'(i ~^wo ) *-®s 4*0 


0 < 4> ^ 4>o 


(77) 


where x^q = ^winglet^^'wing c/Cav is symmetric with respect to 4 >. 
Inserting equation (77) into equation (76) the cp coefficients are 


Cp = — ocld p only 

P Trp Xwo+n-^wo) cos 4 >q ’ 


(78) 


From equation (35) the If,n* coefficients are given by 


n+n*-l 


^nn* “ 




(79) 


p= |n-n*l -1 
odd 
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The distribution of equation (28), for the planar-wing winglet is 
gi ven by 


O.Q = a+e 


<j>o ^ ^ 


TT 

2 


= aCOS y + e + Tpi 0^ 


( 80 ) 


where a is independent of however, e and can be functions of (p. 

Assuming only additional loading conditions, then twist e is zero, and assume 
that the toe in angle ip is constant along the winglet. Then the c^n 
coefficient of equation (35) with equations (77) and (80), becomes 


_ 4 „ >^wo (cos Y + |) 
Can = 


[1 - Xwo(cos Y 


a 


X^o ^ X^o ) cos (^0 


odd n only (81) 


which is the c^n for an untwisted rectangular wing with rectangular winglet 
of chord with winglet extending from span station 

(() = 0 to (|)Q, at constant dihedral cant angle y from the wing plane, with 
toein angle ip. 


Nonplanar Kpn* Coefficients 

For the planar-wing winglet, these coefficients are obtained from 
equation (41). The coefficients in this equation are independent of 
chord distribution, hence, once known for a given <j>o and y, they can be 
combined with any chord distribution for evaluating coefficients. 

Example values of for ({>o = 5tt/32 radians and y = 90 and 75 degrees are 
shown in table I. These values are based on computations with M = 15 and with 
(|)w = wit/16. For an M = 15 computation these Ln*w can be inserted into 
equation (41) and ^nn* computed with (c/cav) given by equations (75) or (77) 
in which, a = cos = cos (wtt/16) and ctq = cos (5ir/32) = .88192. Example 
values of Knn* are presented in table III for a rectangular wing with 
rectangular winglet with = c^.^gi that is, x„o = 1. 
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Loading Solution for Rectangular Wing with = 1 Winglet 
This wing-winglet has the geometry 


Sp “2c Sp ) A 


= _C_ = 

c ’ cav 


1 


( 82 ) 


then, from equation (78) 

'n = Pn* " odd p, n, n* only 


(83) 


Equation (79) becomes 


^ I 

4 ^nn* 


n+n*-.-! 

, 

|n-n*| -1 
odd 


(84) 


and equation (81) reduces to 


Can = ^ tcos Y + J + (1 


cos Y - „)cos n<|)o], odd n only 


Equations (38) and (39) become 

4[C0S Y + ^ + (1-COSy - ^-)C0S (j)o] 


ai - 


A + 2k[l + |Kii+ ^ (l‘'^*l^in*)|^] 

n*=3 

odd 

N 


|a = £aa |!i(l + |k„,)+ f. y[|oa(,+n*|Kin*)-n*f Inn" 


n*=3 

odd 


. n* f K.l |nl 

2 nn*J ai 


(85) 


(86) 


(87) 


for n * 3, 5, 7, • • • N 
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where 


, cos Y + J + (1 - cos Y - ^)cos n<Po 
£ml^1 2 °! ( 88 ) 

cos Y + ^ + (1 - COS Y - ~’)cOS <{)q 


and n*(ir/4) Ipp* values determined from equation (84) are tabulated in table 
IV. 


Slender Wing Theory 

As aspect ratio becomes small, the planform becomes slender and independ- 
ent of wing or winglet planform shape or sweep if trailing edge cutout is 
behind leading edge. With the e = 0, ac distribution of equation (80) 
inserted into the Can integral of equation (45), the Ccm for A 0 is 


^ca ~ - cos Y - J)(i>o - sin ())o cos (j)o)] 

Ce,„ - ^ (1 - COS V -■ . SM^]. odd n > 1 J 

then the ratio is 

(1 - cos Y - [ Si'i(»-^)4>o _ sin(n|l).j)0] 


Can = 


-ai 


J - (1 - cos Y - f) (i-o - sin (|)o cos i>o) 


, odd n > 1 


(89) 


(90) 


Comparing the equation (45) Knn* coefficients for slender wings with the 
Dnn* drag parameter coefficients of equation (67) it can be seen that 


K 


nn 


* 

A -> 0 



(91) 


With equation (91), the slender wing equations (49) for determining symmetric 
loading Fourier coefficients become 


ai - 


‘’^^Cfn ] 



If) 


(92) 
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(93) 


®IL = i ^otn + 1 (Con. 4 p 

ai n n 'cctj ir 



M 

n*=3 

odd 


n*^ D 

Cai 


in*-n*:jj^)nn*) 


an* 

ai 


for n = 3, 5, . . . N 

where Cap is from equation (45), with ac from equation (28), and Dpp* from 
equation (67). Equations (92) and (93) are general for obtaining nonplanar 
slender wing solutions. The first term in equation (93) (that is, with 
Dpn* = 0) is the slender wing planar wing solution derived in reference 8. 

For the planar-wing winglet configuration, the Cdp values are those in 
equations (89) and '(90), and Dpn* coefficients are determined from equation (71). 

an Coefficients for Aspect Ratio Approaching Infinity 

These ap coefficients are given in equation (47) and are directly in 
terms of c^jp. Thus as A -> «> 


ai 



in ^ 
ai 


(94) 


where withcap values of equation (81), the ap values for a rectangular wing 
with rectangular winglet for A are 

^wo (cos y + J)+ [1 - Xwo(cos Y + f)]cos <j)0 
^ A ^wo ( ^ “ ^wo^cos 

> 




-V 


1 t [XiJ(C0S Y + cos n(|)o 

1 + C^wo(cos Y + J)"^-l] cos 4 >q 


(95) 


With Xy/o = 1, equation (95) simplifies tO) 
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ai ^ ^ [cos ((>0 + (1 - cos 4,0) (cos y + ^)] 


( 96 ) 


> 

^ 1 1 + [(cos Y + j^)~^-1]cos n 4»n 

ai 1 + [(cos Y + ^)'-l]cos 4>o ^ 

which is for a rectangular wing with rectangular winglet having the samei 
chord as the wing. 

Rectangular Wing Winglet Numerical Solutions for an 
and Force and Moment Coefficients 

These example solutions are all for rectangular wing and rectangular 
winglet with winglet chord equal to wing chord. However, as aspect ratio 
becomes small these data apply to arbitrary wing and winglet shape 
as born out by the slender wing theory. Parameter values that remain the 
same in this numerical example include 


^0 ~ ^w ~ ^wo ~ ^ ’ '}*o ~ 5tt/32 rad, oq ~ cos (|)q ~ .88192, M - 15, N - 9 


Parameters which have various values include 


(97) 


A, y , i); 


(98) 


With N = 9 equation (87) becomes four linear simultaneous equations in an/ai. 
From equation (24) k = (A + 3.79)/(A + 1.895)- Also define as a shortened 
notation for equation (88) which with y = 90 degrees and 4)0 = 5 /32 radians 
= 28.125 degrees results in 


i + (1 - ^)cos (28.125 n) 

4, = = i « 1: 1 1 

^“1 " .88192 + .11808 ^ 

^ a 


(99) 


Then with I^n* values from table IV and the y = 90 degrees values of K^n* 
from table III the four linear equations become 
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(^+ 2.0702)11^3 - .69737 = (^ + 8.5589 - 1.8682 ,1,3) |^ + (6.5716 - 1.9083 ,^3) 

+ (6.3840 - 2.0005 ,1,3) |^ + (6.2065 - 1.8933 

.44574 = (3.4571 - 1.8682 1/15) + (f + 17.6827 - 1.9083ii-5) 

+ (12.7130 - 2.0005'I»5) (11.5070 - 1.8933 ^ 5 ) 7^ 

a 1 a 1 

( 100 ) 

.40765 = (2.6664 - 1.8682 ^ 7 ) + (9.2030 - 1.9088 ij^y) 7^ 

a 1 a 1 

+ (J-+ 28.0203 - 2.0005 ^, 7 ) 7^+ (18.6383 - 1.8933i|)7) 

(^+ 2.0702)iJ; 9 - .30123 = (2.3797 - 1.8682 ijig) |^ + (6.8409 - 1.9088 ijig) 

+ (14.5977 - 2.0005 <l»9) ~ + (f + 37.2048 - 1.8933iji9) 

a 1 K a 1 

When A and ^/a are specified, a simultaneous solution of these four 

gives the ap/ai from n = 3 through 9. With ap/ai known, aj/a can be evaluated 

from equation (86) which, with y * 90 degrees and <j>o * 5fr/32, becomes 

3.5277 (1 + .13389 ^) 

iL = 2 (lOI) 

“ A + 2k (1.0351 + .9341 + .9542 + .9467 ^) 

ai ai ai a/ 

Loading equations for other y angles using equations (86), (87), and (88) 
can be written in a similar manner to those of equations (100) and (101). It 
can be noted that with y * 0, Kpp* * 0, which is a plannar wing-wing let 
configuration. Values of ap determined from equations (100) and (101) and 
for other Y angles are presented in table V. 

S lender wing . - With A 0, solutions of equations (93) and (92) are 
needed. For if.Q = 28.125 degrees, equation (90) becomes 


(^ + 2 . 0702 ) 1/^5 - 
(^ + 2.0702)'I^7 - 
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Pncn. = 


c 


ai 


- .6686 (1 - cos Y - i ~sin (n-1 ) 28.125 
1 + .05024 (cos Y + ^) ^ 

a 


sin (n+1 ) 28.125j (io2) 


where the value within the brackets for n = 3, 5, 7, . . . odd is respectively, 
.18477, .19846, .12090, .00969, -.09870, -.07157, -.02282, .03164, and so on. 
The coefficient C(jj is determined from equation (89) which, for (j)o = 28.125 
becomes 


1.212334 [1 + .05024 (cos y + ?)] (103) 

a 'a 

then, from equation (92), the ai/a coefficient for A 0, (fro = 5ir/32 radians, 
and N * 9, is 


a 


.95217 [1 + .05024 (cos y 

5 ~ 



TT 


1 1 




ai 


(104) 


For an A ^ 0, N = 9 solution for a^'s, equation (93) is solved. The value of 
^dcn/^il^ given in equation (102) and n*(4/ir)Dnn* values are in table II. 

With Y = 90® equation (93) is 

-.00111 + .04255 (fro = (-.94063 - .00027 ipo) |^ + (.05642 + .00007 to) 

+ (.01882 + .00051 ifro) (.00965 + .00255 to) ^ 

di 

-.00052 + .02742 ifr^ = (.03267 - .00018 ifro) |^ + (-.98172 + .00004 ifrg) 

+ (-.00323 + .00033 tfrg) |^ + (-.02269 + .00033 ifro) 

.00519 + .01193 (fro = (-.00349 - .00008 (fro) + (-.01700 + .00002 (fro) 

+ (-1.02633 + .00014 (frp) |^ + (-.01645 + .00072 (fro) 

.00351 + .00074 (fro = -.01391 - .01908 |^ + (-.00932 + .00001 

+ (-.99630 + .00004 (fro) 
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where 


'Po "" 


1 - ^ 

ot 

1 + .05024 ^ 
a 


With Y = 90 degrees, equation (104) reduces to 


= 

a 


.92144 (1 + .05024 

^ g 

1 + .00645 ^ - .00162 - .01203 - .05992 

ai ai ai a^ 


(106) 


(107) 


When ip/a is specified, a simultaneous solution of the equations of 
equation (105) evaluates an/ai from n = 3 through 9. The first Fourier 
loading coefficient is given in equation (107). A ->• 0 loading equations for 
other Y angles are obtained similarly to the y * 90 degrees example. Values 
of a^, determined for A -> 0 from equations (105) and (107), and for other y 
angles are presented in table V. 

With slender wing nonplanar theory a reasonable approximation for a^/a 
as seen by examining equation (93) is, since Dnn* values are not large 


^n ^ X. 
ai n Cc(i 


(108) 


where Cc^f^/ai for A -> 0 are given in equation (102). A better approximation 
is to use equation (108) for approximating an*/ai under the summation sign. 
Then 

1 *) 


equation U'-’o; ror approxiiiicicing dp 

^ (1 + ^ Dll) - - Dm + - ^ (^ 
Cll T7 IT ni 'Cm 


£an. 
c 


Dm* “ Dnp,7 


roll! 

-oil 


^n /- 
ai 


odd 


n (1 + D 

IT Cctl 


in 


- 7 “nn) 


(109) 


where the prime on the summation sign means the term is not summed when 
n* = n. The coefficients Dip and Dpp symbolize Dm* and Dpp* for n* = n. 

In this form ap/ai can be determined directly and does not involve a 
simultaneous equations solution. When ap/ai are known, ai is computed from 
the equation (92) relation. 
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A •> <” solution . - With <|)o = 28.125 degrees equation (96) becomes 


il. ^ 3.527 7 p ^ ^33p ^ 

an , 1.1339 cos Y + jir + (1 - cos y - cos (28.125 n) 

^ " 1 + .1339 (cos Y + J) 

Equations (110) and (111) are the a^ coefficients for a rectangular wing winglet 
with winglet chord equal to wing chord, as aspect ratio approaches 
infinity. Example numerical values are included in the data of table V. 

Lift coefficient . - With (|)q = 28.125 degrees, N = 9, equation (55) 
becomes 

Cl = n - (1 -cosy) (.047835 + .11763 .12634 .07697 ^ 

2a ' aj ai 3i 

+ .00617 f^)] (112) 

“1 

Using equation (112) with the a^ in table V, the lift-curve slope of the 
various wing-winglet configurations can be evaluated, and these are presented 
in table V. 

Root bending moment coefficient . - With (|)o = 28.125 degrees, N = 9, 
equation (59) becomes 


( 110 ) 

( 111 ) 


^”^^r _ A ai r, + 

a 6a *■' 5 ai 7 ai 15 ai 77 ai 

(113) 

-(1 -cosy) (.099400 + .24443 + .26254 .15994 .01282 ^)] 

ai ai'-" 

Using equation (113) with the an in table V, the wing root bending moment 
coefficient of the various wing-winglet configurations can be evaluated and 
these are presented in table V. 
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Induced drag coefficient . - With N = 9, equation (69) becomes 


n=l n=l n*=l 

odd odd odd 


(114) 


where, with (j)o ■ 28.125 degrees and y angles of 90 and 75 degrees, n*(4/ir) Dnn* 
are listed in table II. Then, with y = 90 degrees, equation (114) leads 
to 


Q 

4^ = 1(^)^ fl. 03335 + (.00333 + 2.82191 ^ - .33261 - .03205 ^ 

+ .09628 ia.) + (-.00425 + 4.90862 + .19985 + .28517 f^) (115) 

dldi ai 

+ ^ (.02387 + 7.18434 + .19903 + fa. (-.03036 + 8.96666 |^)1 

di dii di di Ql J 


An induced drag expression similar to that of equation (115) can be made 
for the Y = 75 degrees configuration by using the D^n* coefficients of table 
II with Y = 75°. Induced drag coefficients computed from equations (114) 
and (115) are presented in table V. 

The drag efficiency factor e is 

_ Cl ^ 

® ''tACq^- COj/irAa^ 

This drag factor e is given in table V as determined from the lift and drag 
coefficient parameters listed in that table. 


Ratio of Induced Drag of Planar-Wing Winglet to that of Flat 
Planar Wing with Equal Lift and Root Bending Moment 

A measure of the efficiency of a wing-winglet configuration is to 
compare the configuration with a flat planar wing with the same lift and root 
bending moment as that of the wing-winglet. The efficiency is measured in 
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terms of induced drag. Let subscript p denote a flat planar wing. Then 
the induced drag is 

Ln^ _ i2 

^ip " irqbp^ep ’ ~ 4irqse^e 

where the second equation is the induced drag of the wing-winglet configura- 
tion. In ratio form 



The ratio of root bending moment coefficients gives 

r u /r. Mbp/2LSg 

Cmbr/Cl Mby. Lp bp 

(Cjnbir/CL )p ^h)"p/LI^p ~ ^b)f’p ^ ^^e 


(117) 


(118) 


For the flat planar wing 


(^)p= 7 ^ ' 1 "CP 

where ncp is the spanwise center of pressure along the semi span. For the 
condition of equal lift and equal bending moment, bp/2se, according to 

equation (118), must be 



and the drag ratio from equation 



The ratio of aspect ratio of the 
is 


Cmbj^ 


( 120 ) 


(118) becomes 


( 


4 


Ticp Cl ^ e (121) 

flat planar wing to that of the wing-winglet 
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(122) 



When surface areas are the same, friction drag will remain about equal. 
With the condition Sp = Sg 




(123) 


then the drag ratio is 



A e 


(124) 


F lat planar elliptic wing . - Since the elliptic wing is the most 
efficient in terms of drag, it serves as a good standard for comparison with 
nonplanar wings. Also the elliptic wing is analytically the simplest. For 
the elliptic wing 


^ 4_ 

'^cp 3 tt ’ 


Sp - 


1 


(125) 


then from equations (123) and (124) 



(3tt 



^ ^ 1 
Dip - r e 


Values of these ratios have been computed using the data of table V. 
Results are shown in table VI. 


(126) 


Flat planar rectangular wing . - Another standard for comparison is to 
compare the nonplanar wing with a flat planar wing having the same chord 
distribution. In the present problem, the nonplanar wing-winglet has a 
rectangular shape. For flat planar rectangular wings, it is convenient to use 


®P 1 + 6 


(127) 
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P, _ 2irAp 

^Lap “ Ap + 2kp(l + t) 


028) 


where, from equation (286) of reference 6 

6 = .01453 lOQe [l + (.52 ^)2 + (.26848 ^)6 + (.07 (129) 

T = .10467 loQe [l + .385 (.19167 (130) 


where the subscript p indicates that the A and k are for the flat planar 
rectangular wing. The value of kp» (Ap + 3.79)/(Ap + 1.895) is given in equation 
(24). Based on the analytical data of the rectangular wing (ref. 6), the 
spanwise center of pressure can be formulated as 


’icp 


1 _ 1.5573 

2 (^+A + 20.6)p 


In a squared formulation 


n 


2 - 

cp -A 


25 


An + 

+ 12.1 


2,183 


21 


(131) 


(132) 


Inserting this ncp 
Thus 


in equation (123) leads to a quadratic equation in Ap. 


Ark ^mby. 

16 


Ap + 12.121 
.25 Ap V 2.183 


hence 


Ar 


[s. 


732 



- 775.744 ( 



(133) 


Thus, with a computed nonplanar wing value of Cfnbp/^L ^ given aspect 
ratio A, the equivalent planar wing Ap is obtained by solving the quadratic 
equation (133). With Ap known, ep is determined from equations (127) and 
(129). Then the drag ratio is determined from equation (124). Values of 
these drag ratios have been computed using the data of table V with results 
shown in table VI. 
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A - As aspect ratio approaches infinity the wing circulation 
becomes directly 


_ o 

^ - 2Trac 


(134) 


then similarly to the formulation of equation (54) and (56) 

1 


/ CXp P 

— - — cos Y da 
“ Cav 


(135) 


Cmbr /■ ^ cos y + 5 sin y) da 

= TT / a Cgy 


(136) 


For the planar-wing winglet, these equations become 


Cl = 2 tt f r da + (cos Y + ^)cOS Y / 7 ^ dal 

*-« y Cav “7 ^av ■' 


(137) 


^br 


= ^ /” 7 ^ crda + ^ (cos Y + — ) f (o - 

1 J cav 2 ' o.' J Cg^ 


oq + Oq cos Y)da 

(138) 


Then for the A ^ «> rectangular wing with winglet chord equal to wing chord 


Cl^ = 2tt [cos 4>o + (1 - cos <{)o)(cos y + f)cos y3 


(139) 


^mb 

— ^ = J [cos2 (|)o + (1 - cos ({)o)(l - cos (|)q + 2 cos <Po cos y) 

(cos Y + ?■)] 


(140) 
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The ratio of equation (140) to (139) gives the bending moment to lift ratio 
as 

'^L 8 [l * * a>'=°^ 

The drag ratio given in equation (124) is the product of the aspect 
ratio and the drag efficiency factor ratio. As A the solution of equation » 

(133) gives 


where Cmb^CL is given in equation (141). For a planar rectangular wing 

1 


(142) 


en - 


y I 

n^i " 

odd 


(143) 


For the nonplanar wing the drag efficiency factor is 


e 


Cl 


a 


ttA Cd^- 


where Cp^./a^ is obtained from equation (114). The A -)■ «> an coefficients are 
presented in equation (96). Combining these equations results in 


1 + 1-cos «{»n 
cos ^0 

Op. = ( ^ 

e ' ^ 1-cos (|)n 

cos (^0 


)2 

ipyCos y 




1 


odd 


1 

n 


( 144 ) 
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where ipy is shortened notation for 


= cos Y + f (145) 

The summation in the denominator of equation (144) is a logarithmic infinity. 

The Dnn* summation in ratio to this is negligibly small. Using the a^ 

coefficients of equation (96), the first summation in the numerator becomes 


W UJ 


n=l 

odd 


n=l 

odd 


^ [’I'y + (1 - ^ ^ 


n=l 

odd 


n=l 

odd 


Tip + (1 - L)cos <PoJ'^ 


Z cos n (pn 
n 


(1-lpy) 


n=l 

odd 


” 1 
I ^ 


n=l 

odd 


cos 2n j)o 

Z_ n 


+ ^ (1 - J (1 - Ipy) 


n=l 
2 odd 
Y' 00 


1 1 


n=l 

odd 


Now 


CO 

^ cosj_ia , I In loot |a| , and ^ ^ |c„t 


i>ol 


n=l 

odd 


n=l 

odd 
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which are finite values for 4 >q > 0. Since these terms are divided by 
infinity the expressions with the cosine summations are hence zero. Equation 
(144) now simplifies to 


ep _ 1 - 2ii)y + 

' ' 2 *0 n - %cos 


(146) 


With equation (141) inserted into equation (142), the planar rectangular wing 
aspect-ratio ratio as A-*-co becomes 



Cmb 


( 8 -q ^)2 = C0s2 (),o' 


L' COS (l>o 


T _L 1 cos ({)fj , 

’ ' COsl,o '^T 


(147) 


Combining equations 
is 



- (1-2 ipY+ 


(124), (146), and (147), the induced drag ratio parameter 

2 






(148) 


where ipy = cos y + ip/a. Equation (148) is the ratio of the induced drag of a 
rectangular planar-wing winglet as A - 4 - «> to the induced drag of a planar 
rectangular wing of equal root bending moment. 


In equation (148), there is a value of \p/a at which Di/Dip is minimum. 
Take the derivative of D-j/Dip with respect to i|y, set to zero, and solve! 

for illy, to obtain 


’^'Ymin 

Di 


^) 

a ' 


mm 

Di 





mm 

Di 


2(1 - T2) 

+ 4(1 - T2)(3 + 3T2 + 6 Tcosy 

cos Y 



> 


(149) 
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where 


T = 1 - cos (frn 
COS 4>o 

p = 3 _ 3J2 _ (3T + t3) cos Y - 2J2 COS^y 
The example configuration has <|)o = 28.125 degrees. Then 


( 150 ) 


T = .13389 

P = 2.94622 - .40406 cos y - .03585 cos^y 


> 


lb » 

' a>niin- ' 

t>1 t>1 

and equations (148) and (147) become 


1.96415 


P + (P2 - 3.92830P + 11.99614 + 3.15571 cos y) 


1/2 


(151) 


(152) 



(1 - + 3lljy^) 


[1 + (.01793 + .26778 cos y)^;^]^ 
(1 + .13389 Tj/y cos y)‘* 



,77779 


1 + (.01793 + .26778 cos y 


1 + .13389 cos Y 




(153) 

(154) 


Using equations (151) through (154), numerical values can be determined. These 
values are presented in table VII. 

An important wing-winglet condition is to have minimum drag ratio occur 
when tp/a is zero. This means that the drag will remain at the bottom of the 
drag bucket at all values of a; that is, independent of lift. Minimum drag 
results when the conditions of equation (149) are fulfilled. When ip/a = 0, 
’^Tmin Di t» then solving equation (149) for cos y leads to a fourth 

degree polynomial in cos y, given by 

2T2cos‘*y + (9T-2T2 +T3) cos^y - (3T-6T2 +T^) cos^y + 3(1-T2 )cos y - 1+T2= 0 

(155) 
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For the example configuration <j>o = 28.125 degrees, then T = .13389, 
and the solution of equation (155) for cos y results in 


cos Y = .32987 


Y = 70.74 degrees 



.80951 


> (156) 


in which Di/Dip and Ap/A are evaluated with ipy = cos y in equations (153) 
and (154). 

Slender wing with y = 0 . - Since D^n* coefficients are zero when y = 0» 
then from equations (90), ,(92), and (93) for the slender wing 


M = 1 
ai n 


a 


1 + 


2. rsin(n-l) 4>n _ sin(n+l) (|)o1 ± 
ir . n-1 n+1 . g 

1 + f (i'o - sin (|)o cos <Po) ^ 


^ ((J .0 - sin (j)o cos (|)o) ^ 


(157) 


(158) 


For Y = 0, equations (55), (59), and (66) with equation (157) can be written as 


Cmbr _ 1 / 

C|_ ~ 3tt I' 



e 


1 + 


1 



g 

1 + ko 



a 


(159) 


(160) 
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where 


ko = I (<l>o -sin <l>o cos (|)o) 
n+1 

(-1 ) 2 rsin(n-l) <j)n sin(n+1 ) <i>o 
(nM) n L n-1 ' n+1 

n=3 
odd 



(161) 


. _ / 2 x 2 1 rsin(n-l) (|)n sin (n+1) <|)o 

^2 ■ n L n-1 " n+1 

n=3 
odd 


As aspect ratio approaches zero, ricp becomes 4/3n, and ep = 1 , then equations 
(123) and (124) for A 0, y = 0, but for arbitrary values of (f>o and \f,/a, 
become 



(37T 



(1 + 


3ki ^ 

^ g 

1 + ko 


2 
ife. 

a 



(1 


+ 


3^' I 
1 + 



0^ 

1 + ko 



(162) 


(163) 


This function in equation (163) results in a drag bucket along ^/a. Take the 
derivative of Di/Dip with respect to ip/a^ set to zero, and solve for \p/a, 
gives 


^) • 
a^min 


1 


i^r 

ki L 


-1 

1 + (1 - 72 


1^ 

k2 


.)V2j 


(164) 


+ kf 
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Inserting equation (164) into (163) gives the equation for minimum drag as 




k, 18 
180 ki'^" 5 


Kil 

k2 


+ 


(2 , h 

^ 180 ki^ 


•)(1 - 72 



( 165 ) 


where ki and k 2 are functions of only (f»o, as can be seen in equation (161). 
With ({>0 = 28.125 degrees, the ko, kj, ^2 from equation (161) are 

ko = .047835, kj = .006857, k2 = .009080 

then from equations (164), (165), and (162) 


a 'min 
Di 


-2.25573 



= .95084 




Cmbr 

= .89868, = .10058 

Cl 


(166) 


For the condition y = 0, ((>0 = 28.125 degrees, values from equations 
(157), (158), (159), (160), (162), and (163) are presented in tables V and 
VI. 
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RESULTS AND DISCUSSION 


Load-Line and Slender Wing Theory 

Analytically, these are two-dimensional theories. That is, only induced 
velocities normal to the surface in the yz-plane enter into the solution. 
Longitudinal effects enter the problem only insofar as wing chord distribution 
in the load-line method influences the solution. The load-line solution is 
for wings with the one-half chord line straight or unswept in the yz-plane. 
However, swept wing parameters, developed in a following section, approximate 
swept planform geometry effects for solution with load-line theory. Because 
of the quasi two-demens ional nature of these theories, they are ideally 
suited for the more complex type of loading solutions (such as, taking into 
account nonplanar effects, titerference effects of multiplanes, fuselage-by- 
image methods, wing-in-jet, ground, tunnel walls, or passing through trailing 
vortices of another aircraft). This is because the load-line method is quasi 
two-dimensional, hence, in many cases, the additional normal -to-surface induced 
velocity due to the added configuration likewise may be determined for two- 
dimensional flow. Because the boundary conditions are satisfied at all points 
with the load-line and slender wing theories, the spanwise loading can be 
predicted accurately for any complex spanwise distribution of wing surface 
inclination (such as winglet toe-in angles, control surface deflections, twist, 
induced angles due to other wings, fuselage, wing-in-jet, ground, or free 
vortex). The introduction of the induced-angle parameter k of equation (24) 
transforms lifting line theory effectively into a lifting surface theory, or 
load-line theory, for the prediction of spanwise loading characteristics. The 
integrated lift is accurate since the definition of k depended on this lift. 

Swept wi ng approximation for load-line theory . - Load-line theory can be 
made applicable to swept wings by a redefinition of the aspect ratio and wing 
chord distribution parameters. From simple sweep theory concept of flow 
normal to wing, the effective aerodynamic semi span of the swept wing is 
given by the distance along the midchord line from wing root to wing tip 
denoted by Sey^. Then the effective aspect ratio and lift parameters are 
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(167) 


A 


A 





Cl 


Equation (167) applies to a sweep angle that varies along the span, such as that 
for a cranked wing. For a constant sweep angle, Sqj^/sq * cos A, then 

= A/cos A; CLy^ = Cl/cos a (168) 


where A is the sweep angle of the wing one-half chord line. 

In addition to the swept-wing parameters of equations (167) and (168), the 
swept wing needs a spanwise wing chord factor which relates wing sweep angle 
to wing taper ratio. Shown in figure 21 of reference 9 is a curve of taper 
ratio as a function of sweep angle. Wings with this sweep and taper ratio 
have approximately elliptical spanwise loading independent of aspect ratio. 

This curve was analyzed for the present study with the results 


X = 


0 - 7 C 1 - sin A 
1 + sin A 


(169) 


Equation (169) forms the basis for defining the equivalent chord distribution 
of a straight wing which will give a spanwise loading distribution similar to 
that of the swept wing. The effective chord distribution is of the form 


/ c \ _ 1 - sin A + 2 |q| sin A c 

'^av e ^ ^ ^ ^ *^cc ^ ^av 


(170) 


where the center of spanwise chord distribution is 

1 


^cc 


/ 


Cav 


ellipse yj 

Stt 


da 


straight 

1 + 2X 
3 + 3a 


(171) 
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When a = 1, c - ct. and, when cj « o, c = c^. These values are inserted into 
equation (170), and a ratio made; then equation (170) becomes equation (169) 
when Xe = .375. Equations (167), (168), and (170) provide swept wing load- 
line theory replacement parameters which are read as 


for A in load-line theory, substitute Sqa A/Sg or A/cos A 

for C|_ in load-line theory, substitute or Cl/cos a > (172) 

for c/cav in load-line theory, substitute {c/Cg^y)^ of equation (170)J 


The chord conditions given in equation (170) provide an interesting 

relationship for swept wings which have elliptical additional spanwise loading 
distribution. For the unswept elliptic planform, (c/cav)g = (4 /tt)( 1 
Then the chord distribution for swept wings is 


e _ 7 [1 - (1 - jf) s1n A] 

Cav ~ 1-sinA +2|:rrsinA 


(173) 


For the chord distribution given in equation (173), the swept -wing aerodyanmic 
characteristics are given by the elliptic wing results of equation (46), 
provided A is replaced by A/cos A and Ci_ by Cl/cos a. 


Rectangular Wing Winglet, (pQ = 5ir/32, x^q = 1 

Example solutions for the aerodynamic characteristics of this configura- 
tion are shown in the previous chapter. Results, including influence 
coefficients and force and moment coefficients, are presented in tables I 
through VII. Examination of table VI shows that, compared to a planar elliptic 

wing of equal lift and equal wing root bending moment, the drag ratio for the 
wing-winglet is about 5% less when A -»■ 0, about 10% greater for A = 10.94, 
and infinitly greater as A Comparing with a planar wing with same 

rectangular planform shows definite drag buckets at all aspect ratios. The 
induced drag ratio, ratio of aspect ratio, and wing root bending moment data 
of tables VI and VII are presented in figure 2 as functions of winglet toe-in 
angle. The induced drag ratio curves of figure 2(a) have drag buckets with 
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respect to toe-in angle. These buckets are shallow and wide for the low 
aspect ratio wing but become deep and narrow as aspect ratio becomes large. 

At minimum induced drag ratio, the comparative reduction of induced drag is 
about 4% at A 0, 11% at A = 10.94, and 66% at A The position of the 

drag bucket along the axis in figure 2(a) is strongly influenced by the 
winglet dihedral cant angle y. It is also influenced by the ratio of winglet 
chord to wing chord. This characteristic is useful for configuration design 
since, when a value of and winglet chord are specified, y is defined 
by the minimum induced drag condition. 

Curves of ratios of wing root bending moment coefficient to lift 
coefficient are presented in figure 2(b). These values apply to both the wing- 
winglet and planar-wing configurations since the condition for comparison was 
that lift and root bending moment be the same. The Ci^bp/Cl 'f'atio becomes 
smaller as aspect ratio and toe-in angle decrease and as dihedral cant angle 
increases. The aspect ratio of the planar rectangular wing in terms of that 
of the rectangular wing-winglet is given in figure 2(c). This is the aspect 
ratio the planar wing must have in order to give the same lift and wing bending 
moment as the wing-winglet with higher aspect ratio. It can be recalled that 
the wing-winglet aspect ratio definition is based on total spanwise distance 
along wing and winglet and total surface area which includes wing area plus 
winglet area. The values of figure 2(c) correspond through A, and y with 
those of figures 2(a) and (b). Figure 3 results from cross plotting the values 
of figure 2. The decrease in wing root bending moment coefficient as the 
winglet approaches the up position is shown clearly in figure 3. For A = 

10.94 at Y = 90 degrees the root bending moment is 4.3% less, however, induced 
drag ratio is 3.4% more than that at y = 0 degrees. The planar rectangular wing 
induced drag D^-p is not a constant in the Di/Dip ratio but varies with the 
planar wing aspect ratio Ap. The value of Ap was specified under the 
conditions that lift and wing root bending moment of the planar wing equal 
that of the nonplanar wing. Curves of only the rectangular wing winglet 
induced drag versus wing root bending moment are presented in figure 4. The 
induced-drag parameter is e"^ = ttA CDi’/Cl^, aspect ratio is 10.94, and the lift of 
the various configurations is the same. As shown in figure 4, at low values of 


58 



root bending moment, the induced drag is marginally less at high values of 
dihedral cant angle. The flat winglet y = 0 configuration curve almost 
envelopes the higher y angled curves. However, the y = 0 curve is characterized 
by large negative values of i|)/a which become less negative as angle of 
attack increases resulting in larger root bending moment coefficients. 

A deep and narrow drag bucket means that induced drag is sensitive to 
changes in ip/a. Thus, when a design toe-in angle is specified, any change in 
angle of attack or lift will alter the ^/a ratio and induced drag rises up 
the side of the drag bucket. In this regard, the ideal would be zero toe-in 
angle; then the i|)/a ratio remains the same at any lift. The dihedral cant 
angle y has a valuable asset in that it influences the position of the induced 
drag bucket along the iji/a axis. Then y can be chosen such that minimum 
induced drag ratio can be obtained with minimum ratio. Examination of 
figure 2(a) indicates that, for this wing winglet configuration, the minimum 
induced drag for \\>/a = 0 occurs at y = 71 degrees as A oo, and at y = 90 
degrees for A = 10.94. In view of this characteristic a good A = 10.94 
rectangular wing winglet has y = 90 degrees, » 0, and minimum induced 
drag. For example, a designer has proposed an aspect ratio equal to nine 
planar rectangular wing, and asks if drag can be reduced without altering 
wing root bending moment. From figure 2(c), when Ap = 9, iji/a = 0, y = 90, 
then A = 10.94. The bending moment coefficient from figure 2(b) is Cmb^/CL = 
.1041 which leads to the same root bending moment as that of the planar 
rectangular wing. The induced drag ratio from figure 2(a) is Di/Dip = .899. 
Thus, the designer has an alternative configuration which has 10.1% less 
induced drag, the same wing root bending moment, and equal surface area. 

Since surface areas are the same, friction drag is approximately the same 
for planar wing and wing winglet configurations. This wing winglet config- 
uration is with winglet starting at i))© = 5 t 7/32 span station. The effect of 
other winglet sizes needs to be investigated for establishing an optimum. 

The swept wing approximation equations can be applied to the rectangular 
wing winglet example. In the rectangular wing solution, (c/cav)g of equation 
(170) was unity. Then 
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(174) 


^ 1 - sin A + 2 acc sin A 

Cav ” 1 - sin A + 2 |a| sin A 


where acc Is given in equation (171). Inserting equation (174) into equation 
(171) and solving for a^c gives 


<^cc 


1 

2 tanh~^sin A 


1 - sin A 

2 sin A 


(175) 


then equation (174) becomes 


c _ (tanh~^ sin a)~^ sin A 
Cav 1 - sin A + 2 |al sin A 


(176) 


where A is the half chord line sweep angle. As an example with A = 30 
degrees 


c ^ 1 .8205 
Cav 1 "t" 2 a 


(177) 


which is not a straight-tapered wing. However, a mean taper ratio is 
approximately X = 0.29. From equation (172) with A = 30 degrees and 
A = 10.94 the swept wing aspect ratio is A = 9.47, and swept wing lift 
coefficient is 0.866 times that of the rectangular wing winglet. Thus, in 
this example, the A « 10.94 rectangular wing winglet results of figures 2, 
3, and 4 apply equally to a 30 degree sweptback wing winglet of aspect 
ratio 9.47, with approximately a taper ratio of 0.29. 


CONCLUDING REMARKS 


A goal for the nonplanar wing loading theories developed in the 
present report was to satisfy boundary conditions of the wing analytically 
at all points. This goal was attained by Fourier analyses which permits 
integration of wing chord and wing inclination boundary conditions which 
specify no flow through the wing. The nonplanar load line, slender wing. 
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elliptic wing, and infinite aspect ratio limit theories have proved to be 
useful tools for applications in loading aerodynamics. They give lifting 
surface accuracy with lifting line theory simplicity, that is, they remain 
quasi two-dimensional solution methods. Because the boundary conditions are 
satisfied at all points, the spanwise loading circulation can be predicted 
accurately for complex spanwise distribution of angle of attack such as that 
due to winglet or control surface deflections, wing twist, induced angles by 
multiwings, multiwinglets, ground, walls, jet, or fuselage. Because of the 
quasi two-dimensional nature of these theories they are well suited for 
these complex type of loading solutions, since, in many cases, the induced 
angles normal to the surface due to the added configurations are likewise deter- 
mined for two-dimensional flow. Force and moment coefficients of nonplanar wings 
involve integrations along the nonplanar surface of loading circulation vectors 
relative to the surface. Equations for nonplanar wing lift, running wing 
bending moment, rolling, and induced drag coefficients are developed. The 
nonplanar wing induced drag coefficient is equal to the planar wing induced 
drag expression plus a summation of products of wing circulation Fourier 
coefficients and nonplanar induced drag influence coefficients. These 
influence coefficients are functions of the nonplanar curve of the wing but 
are independent of any other wing geometry parameter. Swept wing parameters 
are approximated for wing chord and aspect ratio.. These are used in the load 
line method for obtaining nonplanar swept wing solutions. 

Example application of these theories have been made for a rectangular 
wing with wingtip winglets. In this example configuration the winglet extends 
from span station .882 to 1, and winglet chord equals wing chord. The thirty- 
degree swept wing planform equivalent to the rectangular wing has a taper 
ratio of about 0.29. Compared with a planar rectangular wing of equal lift 
and root bending moment, the winglet wing has minimum induced drag ratios of 
.948, .869, and .336 at aspect ratios of 0, 10.94, and «>, respectively. The 
study of this example shows that there is just a deep a drag bucket with 
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winglet planar, but with negative deflection, as with the winglet up at high 
dihedral cant angle. However, it was found that dihedral cant angle has a 
valuable asset in that for the optimized design configuration the drag at 
the bottom of the drag bucket can be realized for any wing angle of attack 
or lift, whereas for the planar winglet the drag rises up the side of the 
drag bucket as angle of attack increases or decreases. 

Vought Corporation Hampton Technical Center 
Hampton, Virginia 23666 
March 21, 1977 
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TABLE I. - COEFFICIENTS FOR EQUATIONS (41) AND (71) OF 
PLANAR-WING WINGLET WITH <f«o= 5 t:/32 aND M=15 


*v 

deg 


1 

2 

3 

4 

5 

6 

7 

8 

90 

1 

.01931 

.13347 

-.11359 

. 04958 

.03096 

.01885 

.01361 

.01221 


3 

-.04873 

-.16124 

-.10834 

.03594 

. 02439 

.01516 

.01104 

.00993 


5 

.04879 

-.05402 

-.09461 

.01476 

.01390 

.00921 

.00687 

.00624 


7 

.02698 

.07104 

-.06937 

-.00490 

.00349 

.00319 

.00261 

.00246 


9 

-.02714 

.06417 

-.03415 

-.01567 

-.00341 

-.00101 

-.00042 

-.00025 


11 

-.02089 

-.01173 

.00316 

-.01888 

-.00542 

-.00256 

-.00163 

-.00135 

75 

1 

.00066 

.08511 

-.06993 

.03114 

. 02046 

.01291 

.00951 

.00858 


3 

-.02812 

-.11212 

-.06670 

.02232 

.01600 

.01033 

.00769 

.00696 


5 

.03297 

-.02911 

-.05830 

.01264 

.00891 

.00627 

.00473 

.00433 


7 

.01538 

.04714 

-.04294 

-.00384 

.00192 

.00202 

.00173 

.00165 


9 

-.01795 

.03687 

-.02156 

-.01048 

-.00262 

- . 00086 

-.00039 

-.00025 


11 

-.01296 

-.01042 

.00110 

-.00982 

-.00381 

-.00186 

-.00120 

-.00100 


TABLE II. - n*^ Dpn* COEFFICIENTS FOR EQUATIONS (69) AND (93) 
OF PLANAR-WING WINGLET WITH 4>o = 5^/32 


Y 

\r\^ 

1 

3 

5 

7 

9 

11 

deg 


. 






90 

1 

.03335 

.00666 

-.00167 

-.01243 

-.06192 

-.10534 


3 

-.00333 

-.17809 

-.16927 

-.05646 

-.02895 

-.06658 


5 

-.00258 

-.16334 

-.09138 

. 08085 

.11343 

.03307 


7 

.03630 

.02441 

.11900 

.18434 

.11512 

-.00789 


9 

.03156 

.12523 

.17174 

.08391 

-.03334 

-.07625 


11 

-.02442 

.06277 

.04171 

-.11282 

-.16333 

-.02769 

75 

1 

.02254 

.01464 

.00887 

-.00752 

-.04600 

-.06860 


3 

-.02833 

-.11801 

-.09743 

-.03368 

-.02510 

-.04331 


5 

-.02362 

-.10743 

-.05579 

.05572 

.06810 

.0:1018 


7 

.02184 

.01044 

.07264 

.11569 

.06939 

-.01221 


9 

.01822 

.08039 

. 10960 

.04817 

-.02089 

-.03748 


11 

-.01661 

.04548 

.02751 

-.07520 

-.09798 

-.00506 
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TABLE III. - n* ^nn*COEFFICIENTS FOR EQUATIONS (39) AND (87) OF 
RECTANGULAR WING WINGLET WITH = 5 tt/32 AND = 1 


y 

N. n* 

1 

3 

5 

7 

9 

11 


n\ 







90 

1 

.03511 

-.06588 

-.04584 

.00025 

-.05336 

-.12659 


3 

.01535 

-.32053 

-.22541 

.01423 

-.00065 

-.13283 


5 

.02287 

-.30001 

-.09515 

. 20900 

.15731 

-.04113 


7 

.06097 

-.02873 

.21046 

.32421 

.12296 

-.08051 


9 

.03951 

.15523 

.31146 

. 14626 

-.12323 

-.14249 


11 

-.04557 

.11283 

.14343 

-.17091 

-.34193 

-.06549 

75 

1 

.02294 

-.04218 

-.01601 

.00123 

-.04333 

-.08389 


3 

.00757 

-.20220 

-.12487 

.01266 

-.01483 

-.09166 


5 

.01103 

-.19751 

-.05132 

.13694 

.08803 

-.03467 


7 

.03396 

-.02045 

.13428 

.20153 

.06996 

-.05948 


9 

.01968 

1 

. 10354 

. 19886 

.08172 

-.07742 

-.07450 


11 

-.03292 

.08485 

.09031 

-.11774 

-.20657 

-.02172 


TABLE IV. - n*^Inn* COEFFICIENTS FOR EQUATIONS (39) AND (87) 
OF RECTANGULAR WING 


\n* 

n\ 

1 

3 

5 

7 

9 

11 

1 

1 

1 

1 

1 

1 

1 

3 

1/3 

4.6 

3.3809 

3.17778 

3.10390 

3.06838 

5 

1/5 

2.02857 

8.93651 

6.14748 

5.59620 

5.37314 

7 

1/7 

1.36191 

4.39105 

13.68594 

9.19620 

8.22020 

9 

1/9 

1.03463 

3.10900 

7.15260 

18.72562 

12.46581 

11 

1/11 

.81368 

2.44234 

5.23103 

10.19930 

23.98962 



TABLE V. - EXAMPLE LOADING CHARACTERISTICS SOLUTIONS OF 


RETANGULAR WING WINGLET WITH AND = 1 


Y 

1 

¥ 

A 






Cmbr 

CDi 

deg 

a 

A 

9] /a 

93/91 

95/91 

a/Zai 


•-a a 


90 

-1 

0 

.8761 

-.0978 

-.0634 

-.0168 

.00561 

1.339A .1304A 

.2081 .8733 



10.94 

.2480 

-.1378 

-.2178 

-.1059 

.0127 

4.263 .4261 

.0217 .7094 



->• oo 

3.055/A 

-.3508 

-.6666 

-.5449 

-.2299 




-.25 

10.94 

.2640 

-.0168 

-.1124 

-.0569 

.0087 

4.4100 .4546 

.0195 .8440 


0 

-»■ 0 

.9217 

-.0463 

-.0317 

-.0057 

.0011 

1.393A .1379A 

.2219 .8857 



10.94 

.2690 

.0236 

-.0772 

-.0405 

.0074 

4.446 .4626 

.0195 .8608 



00 

3.528/A 

.0370 

-.1753 

-.1550 

-.0366 




1 

-> 0 

.9672 

.0003 

-.0029 

.0043 

-.0031 

1.447A .1453A 

.2418 .8769 



10.94 

.2901 

.1614 

.0431 

.0154 

.0028 

4.615 .4963 

.0235 .7690 



-> oo 

4/A 

1/3 

1/5 

1/7 

1/9 



75 

-1 

0 

.9038 

-.0786 

-.0491 

-.0210 

.0064 

1.387A .1372A 

.2151 .9066 



10.94 

.2547 

-.0889 

-.1754 

-.0927 

.0116 

4.348 .4421 

.0203 .7877 



00 

3. 178/ A 

-.2395 

-.5255 

-.4329 

-.1743 




-.5 

^ 0 

.9261 

-.0548 

-.0342 

-.0155 

.0044 

1.416A .1413A 

.2224 .9136 



10.94 

.2652 

-.0135 

-.1091 

-.0610 

.0097 

4.460 .4639 

.0194 .8674 



oo 

3.414/A 

-.0467 

-.2814 

-.2392 

-.0783 




0 

-> 0 

.9463 

-.0311 

-.0194 

-.0100 

.0025 

1.441A .1451A 

.2294 .9154 



10.94 

.2750 

.0620 

-.0427 

-.0293 

.0078 

4.581 .4905 

.0198 .8964 



^ 00 

3.650/A 

.1211 

-.0689 

-.0705 

.0053 



0 

-2.256 

^ 0 

.8911 

-.1002 

-.0645 

-.0281 

-.0018 

1.400A .1407A 

.2097 .9466 


-2 

10.94 

.2475 

-.1037 

-.1905 

-.1015 

.0100 

4.251 .4320 

.0197 .7772 


-1.25 

10.94 

.2651 

-.0064 

-.1038 

-.0582 

.0089 

4.553 .4863 

.0189 .9272 


-1 

^ 0 

.9522 

-.0412 

-.0265 

-.0116 

-.0007 

1.496A .1553A 

.2288 .9906 



10.94 

.2706 

.0260 

-.0750 

-.0438 

.0086 

4.648 .5045 

.0191 .9577 



oo 

3.528/A 

.0371 

-.1753 

-.1550 

-.0365 




0 

0 

1 

0 

0 

0 

0 

. 5 ttA A /6 

1/4 1 



10.94 

.2914 

.1556 

.0405 

.0139 

.0071 

5.005 .5779 

.0230 .9236 



oo 

4/A 

1/3 

1/5 

1/7 

1/9 
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TABLE VI. - COMPARISON OF INDUCED DRAG OF RECTANGULAR WING 


MINGLET, ((>0 = 5tt/32,. ^wo = ^ * WITH PLANAR WING OF 
EQUAL LIFT AND EQUAL WING ROOT BENDING MOMENT 


Y 

deg 

a 

A 

e 

Cmbr 

“cr 

Planar ellip- 
tic wing, ep=l 

Ap/A Di/Dip 

Compared with planar 
rectangular wing 

ep Ap/A Di/D|p 

90 ' 

-1 

0 

.8733 

.0974 

.8421 

.9642 

1 

.8421 

.9642 



10.94 

.7094 

.1000 

.8874 

1.2509 

.9470 

.7664 

1.0231 


-.25 

10.94 

.8440 

.1031 

.9440 

1.1185 

.9393 

.8111 

.9028 


0 

0 

.8857 

.0990 

.8702 

.9825 

1 

.8702 

.9825 



10.94 

.8608 

.1041 

.9618 

1.1173 

.9380 

.8252 

.8993 


1 

0 

.8769 

.1005 

.8965 

1 . 0224 

1 

.8965 

1.0224 



10.94 

.7691 

.1075 

1.0270 

1.3354 

.9334 

.8768 

1 .0642 

75 

-1 

-V 0 

.9066 

.0989 

.8687 

.9583 

1 

.8687 

.9583 



10.94 

.7877 

.1017 

.9184 

1.1658 

.9412 

.7907 

.9448 


-.5 

0 

.9136 

.0998 

.8846 

.9683 

1 

.8846 

.9683 




.8674 

.1040 

.9612 

1.1082 

.•9381 

.8248 

.8920 


0 

^ 0 

.9154 

.1007 

.9007 

.9840 

1 

.9007 

.9840 



10^94 

.8964 

.1071 

1.0181 

1.1358 

.9340 

.8698 

.9063 

0 

-2.256 

0 

.9466 

.1005 

.8977 

.9483 


.8977 

.9483 


-2 

10.94 

.7772 

.1016 

.9175 

1.1804 


.7900 

.9568 


_-1.25 

10.94 

.9272 

.1068 

1 .0136 

1 .0931 


.8662 

.8729 


-1 

0 

.9906 

.1038 

.9572 

.9664 

1 

.9572 

.9664 



- 10-.94 

- / 9577 „ 

.1086 

1 .0467 

1.0929 

.9321 

.8924 

.8685 




1 

.1061 

1 

1 

1 

1 

1 



10.94 

.9236 

.1155 

1.1842 

1.2822 

.9236 

1 

1 
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TABLE VII. - COMPARISON OF INDUCED DRAG OF A ^ RECTANGULAR 
WING WINGLET, = 5tt/32, A^o = 1 , WITH PLANAR 
RECTANGULAR WING OF EQUAL LIFT AND EQUAL WING 
ROOT BENDING MOMENT 



MINIMUM DRAG AND CONDITIONS FOR MINIMUM DRAG, A - 


Y, deg 

105 

90 

75 

60 

45 

30 

15 

0 

70.74 

T^/a 

.5582 

.3294 

.0709 

-.1697 

-.3761 

-.5343 

-.6338 

-.6677 

0 

Ap/A 

.7691 

.7870 

.8047 

.8210 

.8364 

.8457 

.8524 

.8547 

.8095 

Cmbp/CL 

.1096 

.1109 

.1121 

.1133 

.1143 

.1150 

.1154 

.1156 

.1125 

Di/Dip 

.3373 

.3373 

.3371 

.3368 

.3364 

.3361 

.3359 

.3358 

.3371 
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(a) Induced drag ratio. 

Figure 2. - Characteristics of rectangular wing winglet, <1)0 = 5ir/32, 

Awo = I5 compared with planar rectangular wing of equal lift 
and equal wing root bending moment. Subscript p refers to 
planar wing value 












^mbr 

Cl 


Figure 3. - Minimum induced drag as a function of wing root bending 
moment coefficient from data of figure 2. The angle. y» 
in degrees, is the winglet dihedral cant angle measured 
from the wing plane. 
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Figure 4. - Induced drag parameter e~^ = ttACdi/Cl^ of A = 10.94 

wing-winglet as a function of wing root bending moment 
coefficient, from data of Table V. 
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